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Abstract 

An  appropriate  oxygen  excess  ratio  could  improve  the  net  power  of  a  fuel  cell  system  and  reduce  the  potential  of  oxygen  starvation.  Due 
to  nonlinear  and  time-varying  characteristics,  it  is  a  difficult  task  to  maintain  the  oxygen  excess  ratio  within  an  ideal  range  particularly  during 
transients.  An  adaptive  control  algorithm  is  therefore  proposed  to  dynamically  stabilize  the  oxygen  excess  ratio  around  an  optimal  level.  The 
application  of  the  algorithm  inevitably  involves  the  estimation  of  time-varying  parameters  and  the  pole  assignment  of  the  closed-loop  system.  The 
proposed  adaptive  control  is  evaluated  on  a  test  bench,  and  its  application  to  overcome  the  time-varying  problem  is  demonstrated. 

©  2008  Elsevier  B.V.  All  rights  reserved. 
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1.  Introduction 

Stimulated  by  attempts  to  overcome  energy  shortages  and 
alleviate  pollution  caused  by  emissions  from  energy  consump¬ 
tion,  many  countries  devote  great  effort  to  exploring  and  making 
use  of  new  energy  resources.  Because  proton  exchange  mem¬ 
brane  (PEM)  fuel  cells  boast  many  ingenious  features,  such  as 
high  power  density,  fast  startup,  suitability  for  discontinuous 
operation,  low  noise,  and  environmental  friendliness  [1-4],  they 
are  now  under  intensive  development  by  many  governments, 
manufacturers,  and  researchers. 

The  air-supply  subsystem  is  very  important  in  a  PEM  fuel 
cell  system  [5].  Excess  oxygen  replenishment  into  the  cathode 
will  cause  power  waste,  consequently  leading  to  a  decrease  in 
the  net  power  of  the  fuel  cell  system.  Conversely,  if  fuel  cells 
are  subject  to  oxygen  starvation,  the  potential  of  these  cells  will 
suffer  from  deterioration.  If  the  starvation  is  not  too  serious  and 
the  oxygen  is  enough  to  maintain  the  stack  current,  the  voltage 
of  the  fuel  cells  is  positive.  However,  if  fluxes  of  oxygen  are 
not  enough  to  maintain  the  stack  current,  there  must  be  an  elec¬ 
trolysis  current  that  makes  the  voltage  of  the  starved  cells  drop 
to  negative-named  “cell  reversal”  [6].  In  this  case,  the  protons 
transported  from  the  anode  will  be  oxidized  in  the  cathode  to 
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form  hydrogen  and  the  cell  essentially  acts  as  a  hydrogen  pump 
[7],  The  electrode  half  reactions  and  the  overall  reaction  are  as 
below: 

anode  :  H2(a)  -+  2H+  +  26- 

cathode  :  2H+  +  2e"  -»  H2(c)  (1) 

overall :  H2(a)  -+  H2(c) 

While  the  process  of  air  replenishment  is  directly  related  to 
oxygen  starvation  [8-10],  many  researchers  have  been  focus¬ 
ing  on  the  air-supply  kinetics  in  the  past  decade.  Vahidi  et  al. 
[11]  investigated  oxygen  starvation  in  a  fuel  cell/super-capacitor 
hybrid  system.  The  proposed  model-predicting  controller  can 
avoid  the  oxygen  starvation  phenomenon  despite  load  fluc¬ 
tuations.  Pukrushpan  [12]  designed  a  model-based  air-supply 
controller,  which  could  optimize  the  air-supply  subsystem  in 
principle  but  this  lacks  experimental  support.  Methekar  et  al. 
[13]  presented  a  distributed  PEM  fuel  cell  model  and  studied 
its  dynamics.  Subsequently,  a  proportional  control  strategy  has 
been  attempted  to  reduce  oxygen  starvation.  Feroldi  et  al.  [14] 
proposed  a  coordinative  control  strategy  that  regulates  the  air¬ 
flow  replenished  into  the  system  and  the  exhaust  exiting  the 
system  simultaneously.  As  a  result,  both  the  oxygen  excess 
ratio  and  fuel  cell  voltage  are  well  under  control.  However, 
the  effectiveness  of  the  proposed  strategy  has  not  been  veri¬ 
fied  with  experiments  thus  far.  Rodatz  et  al.  [15]  developed  a 
dynamic  model  of  an  air-supply  subsystem  as  well  as  a  linear 
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Nomenclature 

A  area  (cm2) 

ca  cathode  of  fuel  cell  stack 

e_  minus  electron 

ex  expectation 

/  gas  flow  (mol  s " 1 ) 

F  Faraday  Constant 

gfr  gas  flow  regulator 

H+  proton 

H2  hydrogen 

i  current  (A) 

N  number  of  cells 

N2  nitrogen 

O2  oxygen 

p  pressure  (kPa) 

R  universal  gas  constant  (8.3143  J  Kmol-1) 

u  voltage  (V) 

V  volume  (L) 

t  time  constant  (s) 


quadratic  regulator  for  a  PEM  fuel  cell  system.  According  to 
the  results  presented  in  this  paper,  the  improvement  in  perfor¬ 
mance  achieved  by  this  controller  was  noticeable.  Moreover, 
the  pressure  control  is  successfully  decoupled  from  the  mass 
flow  control,  resulting  in  a  fast  response.  Pukrushpan  et  al. 
[16,17]  presented  a  single-input  single-output  (SISO)  linear  state 
feedback  controller  for  the  air  stoichiometric  ratio.  Although 
a  semi-empirical  fuel  cell  model  was  used,  it  does  provide  an 
insight  into  the  transient  behavior  of  the  air-supply  subsystem. 
Caux  et  al.  [18]  studied  an  approach  to  control  the  air-supply 
subsystem  which  is  particularly  effective  during  transients.  A 
specific  balance  model  is  used  to  maintain  a  constant  pressure 
on  the  cathode  compartment  and  to  follow  a  desired  airflow  rate. 

The  complex,  time-varying  and  nonlinear  dynamics  of  the 
air-supply  subsystem  is  usually  approximated  by  the  principle 
of  mass  conservation,  theories  of  fluid  kinetics,  thermodynam¬ 
ics  and  heat  transfer  under  various  assumptions  and  constraints. 
Popularly  applied  models  all  inaccurately  describe  the  dynam¬ 
ics,  and  the  accuracies  of  the  approximations  vary  with  the  shift 
of  operating  points.  This  problem  has  attracted  much  attention. 
Yang  et  al.  [19]  proposed  a  systematic  method  of  identification 
and  control  of  a  PEM  fuel  cell.  Supported  by  simulation  results, 
an  adaptive  controller  has  proven  to  be  robust  despite  the  vari¬ 
ation  in  parameters  and  the  presence  of  disturbances.  Golbert 
and  Lewin  [20]  developed  a  time-varying  state  space  model, 
analyzed  its  controllability,  and  presented  a  nonlinear  adaptive 
controller.  These  adaptive  strategies  indicate  a  trend  that  future 
attempts  to  develop  nonlinear  models  must  oversee  the  entire 
operating  range  to  optimize  the  stable  and  dynamic  performance 
of  fuel  cell  systems. 

As  shown  in  previous  literature,  air-supply  subsystems 
exhibit  time-varying  and  nonlinear  properties.  This  inevitably 
poses  a  difficulty  in  accurate  and  dynamic  regulation  of  air 
replenishment.  While  most  works  focusing  on  these  problems 


Fig.  1 .  Fuel  cell  system. 


were  presented  with  simulation  results,  this  paper  presents 
experimental  validations  of  an  adaptive  strategy. 

The  remainder  of  this  paper  is  arranged  as  follows.  The  kinet¬ 
ics  of  an  air-supply  subsystem  and  its  time- varying  and  nonlinear 
property  are  introduced  in  Section  2,  followed  by  a  design  of  the 
adaptive  control  strategy  to  achieve  an  optimal  oxygen  excess 
ratio  in  Section  3.  After  discussing  the  experimental  results  in 
Section  4,  conclusions  are  presented  in  Section  5. 

2.  Dynamics  of  air-supply  subsystem 

This  paper  examines  a  membrane-humidifying  PEM  fuel 
cell,  as  shown  in  Fig.  1,  with  an  active  cell  area  of  150  cm2, 
24  Nafionll2  cells  in  series,  and  1  kW  gross  power,  operating 
temperature  of  65  °  C  and  pressure  in  the  stack  of  30  kPa.  Fig.  2 
shows  the  structure  of  our  fuel  cell  system.  There  are  four  major 
control  subsystem  loops  that  adjust  the  air/hydrogen  supply,  the 
humidity  of  reactants,  the  temperature  in  stack,  and  the  power 
output  [21], 

Membrane  humidifying  is  accepted  as  an  excellent  mecha¬ 
nism  for  PEM  fuel  cells  [22,23].  In  this  configuration,  the  thin 
humidifying  membrane  (thickness  of  0.051-0.183  mm)  makes 
it  possible  for  gas  flow  to  have  sufficient  contact  with  the  sat¬ 
urated  vapor,  thus  humidifying  is  efficient  [24],  Therefore,  we 
assume  here  a  perfect  air/hydrogen  flow  humidifier.  In  addition, 
a  cooler  is  applied  to  regulate  the  excess  heat  generated  during 
reactions.  Due  to  the  fast  dynamics  of  the  cooler,  the  tempera¬ 
ture  of  the  stack  can  be  controlled  to  be  within  a  proper  range. 
Moreover,  a  fast  proportional  controller  of  the  fuel  flow  ensures 


Fig.  2.  Scheme  of  fuel  cell  stack. 
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Atmosphere 


Fig.  3.  Scheme  of  air-supply  system. 

the  pressure  in  anode  compartment  tracks  that  in  the  cathode 
compartment  almost  instantaneously.  All  these  controllers  and 
assumptions  are  by  no  means  trivial  but  make  the  concentration 
of  this  paper  on  the  control  of  the  air-supply  subsystem  possible. 

As  shown  in  Fig.  3,  the  air-supply  system  mainly  consists  of 
a  compressor  driven  by  an  electric  motor,  an  air  tank  for  com¬ 
pressed  air  reservation,  an  gas  flow  regulator  to  dynamically 
regulate  the  replenishment  of  air  into  the  cathode,  and  the  vol¬ 
ume  of  the  cathode  manifold  that  houses  the  air-supply  pipes, 
cathode,  cooler  and  humidifier.  Additionally,  the  air  has  to  be 
humidified  to  prevent  drying-out  of  the  membrane  of  the  fuel 
cells. 

2.7.  Modelling  of  the  airflow  channel 

Before  modelling  the  air-supply  subsystem,  several  assump¬ 
tions  and  simplifications  are  introduced,  so  we  can  focus  on  the 
transient  behavior  of  air  replenishment:  (1)  all  gases  obey  the 
ideal  gas  law  due  to  the  low  operating  pressure  (with  respective 
to  the  atmospheric  pressure);  (2)  as  the  volume  of  the  stack  is 
relatively  small  and  the  coolant  can  sufficiently  carry  the  heat 
generated  in  the  process  of  the  reaction,  the  temperature  of  the 
air  inside  the  cathode  equals  to  the  bulk  stack  temperature  that 
is,  in  turn,  equal  to  the  temperature  of  the  coolant  exiting  the 
stack;  (3)  the  properties  of  the  gas  flow  exiting  the  cathode  such 
as  temperature  and  pressure  are  assumed  to  be  the  same  as  those 
inside  the  cathode  and  they  dominate  the  reaction  in  the  catalyst 
layers  in  the  membrane;  (4)  because  the  humidifier  is  effective 
as  explained  previously,  the  gas  in  the  cathode  is  almost  fully 
humidified  and  the  liquid  water  generated  inside  the  cathode  is 
no  longer  vaporized;  (5)  any  extra  water  inside  the  cathode  could 
be  carried  out  with  the  exhaust  gas;  (6)  the  flow  channel  and  gas 
diffusion  layer  are  lumped  into  one  volume,  i.e.  the  spatial  vari¬ 
ations  are  neglected;  and  (7)  the  gas  pressure  in  the  air  tank  can 
be  kept  constant  on  the  assumption  that  air  could  be  replenished 
into  the  air  tank  by  the  compressor  under  the  regulation  of  an 
independent  controller. 

The  mass  continuity  of  the  oxygen  and  nitrogen  inside  the 


cathode 

volume 

and  ideal : 

gas  law  yields: 

dpo2,ca 

RTci 

df 

~  ~Vc7 

-(/o2,gfr  — 

/o2,used  —  /o2,ca,out) 

(2) 

ca 

RTct 

df 

-  ~v^ 

(/N2,gfr  - 

/N2,ca,out) 

(3) 

where  d  is  the  differential  operator;  /o2,used,  /o2,ca,out  and 
/N2,ca,out  are  mole  flow  rate  of  O2  used  in  reaction,  mole  flow 
rates  of  O2  and  N2  in  stack  outlet,  respectively. 


Experiment  results  [25]  showed  the  gas  flow  regulator  is  a 
typical  nonlinear  system,  which,  given  the  input  and  output  sig¬ 
nals,  can  be  approximated  by  linear  models  at  various  operating 


points  as 

Fgfr  ^  =  ~  /air.gff  +  ku>/^’gffMgff 

(4) 

/02,gfr  =  0.21/airgfr 

(5) 

/N2,gfr  =  0-79/air, gfr 

(6) 

where  kwf  —  1  /(22.4  x  60)  mol  (L  s)-1  and  KgT  is  the  gain  of 
the  gas  flow  regulator. 

Eq.  (2)-(6)  consist  of  the  model  of  the  air-supply  subsystem. 

2.2.  Analysis  of  the  property  of  the  gas  flow  regulator 

As  the  gas  flow  regulator  shown  in  Fig.  4  is  a  component  that 
replenishes  air  into  the  stack  cathode,  its  dynamics  are  critical 
to  the  air-supply  subsystem.  The  gas  flow  regulator  is  a  complex 
of  mechanics  and  electronics.  Its  inputs,  outputs  and  parameters 
are  subject  to  the  variations  in  the  inner  and  outer  environments 
such  as  electromagnetism,  temperature,  humidity  and  gas  pres¬ 
sure  at  the  outlet.  What  is  more,  there  is  an  obvious  nonlinear 
property  in  some  operating  regions  especially  at  a  small  com¬ 
mand  voltage  and  this  property  has  distinct  influence  on  the 
dynamics  of  airflow,  which,  from  control  point  of  view,  can  be 
treated  as  a  time- varying  property  at  different  operating  points. 
For  these  reasons,  investigations  were  carried  out  to  identify  the 
disturbance  and  time- varying  property  of  the  gas  flow  regulator. 

Many  tests  [25]  yielded  an  experimental  model  of  the  gas 
flow  regulator  at  a  certain  operating  point  according  to  Eq.  (4): 

Fgfr  “  =  ,/aii\!2fr+ kwf  KgfrUgfrf- V/ajr,gfr+ kwf  Xgfr VWgfr 

(7) 

where  V/ajr  gfr  and  Vugfr  are  disturbances  and  higher-order 
dynamics  of  /air,gfr  and  wgfr,  respectively.  Experimental  data  [25] 
showed  that  V  /air,gfr  and  Vwgfr  both  have  Gaussian  distributions 
with  zero  means  to  some  extent.  Because  they  are  independent 
of  one  another,  their  summation  still  has  a  Gaussian  distribution 
as  shown  in 

e  =  V /aj r  gfr  +  kjj;/AgfrV Wgiv  e  ~  (0,  <Tgfr)  (8) 

Substituting  Eq.  (8)  into  Eq.  (7),  a  simplified  model  of  the 
gas  flow  regulator  is  obtained  as 

Tgfr  d^dT  =  ~  /ak'gfr  +  kwfK  gfrMgfr  +  e  (9) 


Signal  of  airflow 


Inlet 


Gas  flow 


Reference 


_  I  Outlet 


Control  valve 


Fig.  4.  Scheme  of  gas  flow  regulator. 
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Thus  far,  Eq.  (4)  can  be  modified  as  Eq.  (9)  in  consideration 
of  disturbances  and  higher-order  dynamics  of  /air,gfr  and  ugfr. 

To  investigate  the  time- varying  property  of  the  gas  flow  reg¬ 
ulator,  tests  are  carried  out  to  identify  the  parameters  of  the  gas 
flow  regulator  at  various  operating  points. 

For  convenient  analysis  of  the  time-varying  property  by 
parameter  identification,  the  linear  system  model  Eq.  (9)  is  repre¬ 
sented  as  a  discrete-time  equation  of  the  difference  operator,  and 
this  leads  to  an  auto-regressive  moving-average  with  an  auxiliary 
input  model,  the  ARMAX  model  [26]: 

f  Agfriz  *)/air,gfr(k)  =  Bgfriz  '  )Ugfr(k)  +  e(k) 

<  AgfrCz-1)  =  1  +  a\z~x  +  a2z~2  4 - f  arz~r  (10) 

1  Bgfr(z~l)  =  biz~l  +  b2z~2  H - 1-  brz~r 

where  integer  k  denotes  the  time  instant  of  sampling  and  r  is  the 
system  order.  In  a  succinct  form,  the  term  <4gfr(7~l  )./air,gfrW  is 
the  auto-regression  part,  where  z- 1  is  a  backward  shift  operator, 
and  z_/'/air,gfr(&)  represents  the  past  output  /air,gfr(£  —  h),h  = 
1, 2, . . . ,  r;  Bgft(z~l)uam(k)  denotes  the  moving  average  of  the 
past  inputs,  where  z~huam{k)  =  uam(k  —  h),  h  =  1, 2, . . . ,  r; 
and  the  auxiliary  input  e(k)  is  included  to  better  describe  the 
unmodelled  dynamics  of  the  system,  including  higher-order 
dynamics  and  disturbances.  Therefore,  the  corresponding  trans¬ 
fer  function  between  output  and  input  in  the  discrete-time  form 
becomes 


where 

Ibgfr  —  [—  /air.gfrO  —  1)>  ‘  ■  ■  ,  —  /air.gfrO  —  r), 

Ugh (i  -  1), ... ,  Hgfr(i  -  r)]T 

0gfr  =  [a\,  02, . . . ,  ar,  b\,  b2, . . . ,  br]T 
i=  1,2,  ...,L 
L  =  1200 

where  T  is  the  transposition  operator. 

In  addition,  the  characteristic  of  the  noise,  e(k),  can  be 
described  by  its  first  and  second  order  moments,  expressed  as 

E{-t}  =  [E{e(l)}f  £{e( 2)},  •  •  • ,  E{e{L)}]r  (13) 


CovCt]  =  E{-t~tT}  = 


E{e\  1)} 

E{e(lM2)} 

E{e(l)e(L)} 

E{e(2)e(l)} 

E{e\ 2)}  ■  •  • 

E{e(2)e(L)} 

E{e(3)e(l)} 

E{6(3)6(2)} 

£{e(3)e(L)} 

E{e(L)e(l)} 

£{e(L)e(  2)} 

E{e\L)} 

where  E{ }  is  the  expected  value  operator  and  Cov  is  the  corre¬ 
lation  operator. 

The  criterion  function  can  be  adopted  as 
L 

Jgfr(0gfr)  =  'y  ^[/air,gfir(k)  —  ^gfr(k)0gfr]2  (15) 


Ggfr  — 


%r(7) 

Agfr(z) 


(ID 


According  to  the  theory  of  parameter  identification,  all  sys¬ 
tem  modes  can  be  excited  by  a  white  noise,  whose  frequency 
contents  consist  entirely  of  sinusoidal  waves  of  the  same  ampli¬ 
tude  or  strength.  In  the  real  world,  a  signal  with  a  rich  content 
of  frequencies  can  be  generated  by  an  M-series,  and  this  signal 
approximately  satisfies  the  persistent  excitation  requirement  for 
system  identification. 

In  the  process  of  identification,  the  model  of  the  gas  flow  regu¬ 
lator  is  regarded  as  a  black  box,  as  shown  in  Fig.  5.  The  reference 
command,  Ug[t(k),  is  an  M-series  with  mean  value  at  the  operat¬ 
ing  point,  and  /air,gfr(k)  is  the  airflow  response  at  the  excitation 
of  Mgfr(k).  Ggfr(z_1)  mathematically  represents  the  dynamics  of 
the  gas  flow  regulator.  e(k)  stands  for  the  noise  arising  from  the 
inner  and  outer  environments,  and  other  disturbances. 

For  convenience  of  identification,  Eq.  (10)  is  rewritten  as 


/air,gfr(fc)  =  fcgfr(fc)0gfr  +  e(fc) 


(12) 


Fig.  5.  Blackbox  for  parameter  identification. 


Minimizing  JguiOgh)  gives  the  estimation  of  0gfr. 

The  parameters  of  the  gas  flow  regulator  used  in  the  experi¬ 
ment  are:  a  main  time  constant,  about  2  s;  a  frequency  band,  30 
Hz;  an  airflow  range,  0-100  L  min-1  responding  to  the  reference 
command,  0-5  V.  The  characteristic  polynomial  of  the  M-series 
is  adopted  as 

F{s)  =  s10  <g>  s9  <g>  sA  <g>  ^  <g>  1  (16) 

and  the  period  of  the  impulse  is  (l/60)s;  the  test  time  is  20  s 
and  the  datum  length  is  1200.  The  normalized  M-series  used 
in  the  test  process  is  illustrated  in  Fig.  6.  Before  the  test 
begins,  the  compressor  is  run  under  its  controller  to  replen¬ 
ish  air  into  the  air  tank  until  the  pressure  reaches  200  kPa.  In 
the  process  of  the  test,  the  opening  extent  of  the  purge  valve 
(OTPV)  is  divided  into  six  levels,  and  at  each  level,  M-series 
voltages  are  exerted  on  the  gas  flow  regulator  as  the  refer¬ 
ence  command  with  mean  values  of  1  V,  1 .5  V,  2  V,  2.5  V,  3  V 
and  3.5  V.1  For  each  test  state,2  two  tests  are  performed,  and 
the  input  signals  (M-series)  and  output  signals  (airflow  from 
the  gas  flow  regulator)  are  recorded  simultaneously  for  further 
analysis. 

The  test  data  are  analyzed  by  minimizing  the  criterion  func¬ 
tion  Eq.  (15),  and  the  identified  model  is  transformed  into  the 
continuous  model  that  has  the  form  of  Eq.  (9)  with  the  parameters 
of  Xgfr  and  rgfr  shown  in  Tables  1  and  2. 


1  These  voltages  are  typical  operating  points  in  the  gas  flow  regulator. 

2  Test  state  is  defined  as  the  cross-point  of  each  level  of  OTPV  and  each  level 
of  reference  command  voltage. 
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M-series 


Fig.  6.  M-series. 


Note:  The  symbol,  *,  means  the  test  is  not  ran  to  prevent  the 
gas  flow  regulator  from  risking  excess  pressure. 

From  Table  1,  the  reference  commands  with  mean  values 
of  1  and  1.5  V  correspond  to  Kgfr  values  no  larger  than  3 
L  (minV)-1,  which  are  obviously  smaller  than  those  for  the 
other  reference  commands,  which  have  values  larger  than  19  L 
(minV)-1.  Therefore,  before  the  reference  command  reaches 
an  opening  value  of  about  1.9  V,  the  airflow  is  almost  out 
operation.  Similarly,  the  data  in  Table  2  also  shows  the  differ¬ 
ence  between  reference  commands  with  mean  values  of  1  and 
1.5  V  and  other  reference  commands.  The  average  values  of  xg fr 
for  reference  commands  with  mean  values  of  1  and  1.5  V  are 
almost  half  those  of  xgfr  for  other  reference  commands.  This 
result  demonstrates  that  the  main  time  constant,  xam,  behaves 
discontinuously  in  the  whole  range  of  operation.  In  addition, 
Kam  varies  slightly  on  the  condition  that  the  reference  com¬ 
mand  is  larger  than  2  V,  being  within  the  range  of  19.40  to 
20.40 L  (minV)-1.  However,  xam  undergoes  obvious  variation 
with  change  in  the  reference  command  and  OTPV.  To  fur- 


Table  1 

Estimation  of  K gfr 


Reference  command  voltage  o 

n  gas  flow  regulator 

1.0V 

1.5V 

2.0  V 

2.5  V 

3.0  V 

3.5  V 

T.S.  (1.1) 

1.18 

0.89 

19.81 

* 

* 

* 

T.S.  (1.2) 

1.07 

0.88 

19.59 

* 

* 

* 

o 

T.S.  (2.1) 

1.34 

0.70 

19.41 

19.48 

19.55 

* 

T.S.  (2.2) 

1.42 

0.81 

19.46 

19.56 

19.54 

* 

T 

T.S.  (3.1) 

1.51 

0.81 

20.31 

19.65 

19.99 

19.74 

T.S.  (3.2) 

1.08 

0.81 

19.76 

19.62 

19.53 

19.86 

P 

T.S.  (4.1) 

2.55 

0.96 

20.25 

20.11 

19.74 

20.04 

T.S.  (4.2) 

1.52 

1.00 

20.56 

20.11 

20.68 

21.06 

V 

T.S.  (5.1) 

1.40 

0.77 

20.56 

20.44 

21.08 

20.21 

T.S.  (5.2) 

1.39 

0.81 

20.40 

20.88 

20.02 

20.35 

T.S.  (6.1) 

1.29 

0.89 

20.51 

20.58 

20.57 

20.20 

T.S.  (6.2) 

1.33 

0.85 

20.45 

20.50 

19.99 

20.40 

Table  2 
Estimation  of 

Tgfr 

Reference  command  voltage  o 

n  gas  flow  regulator 

1.0V 

1.5V 

2.0  V 

2.5  V 

3.0V 

3.5  V 

T.S.  (1.1) 

0.78 

0.87 

1.51 

* 

* 

* 

T.S.  (1.2) 

0.89 

0.92 

1.39 

* 

* 

* 

O 

T.S.  (2.1) 

0.87 

0.98 

1.83 

1.85 

1.99 

* 

T.S.  (2.2) 

0.80 

0.98 

1.64 

1.80 

1.73 

* 

T 

T.S.  (3.1) 

0.85 

1.09 

1.96 

2.06 

1.86 

2.36 

T.S.  (3.2) 

1.02 

0.99 

1.94 

1.91 

2.07 

2.06 

P 

T.S.  (4.1) 

1.90 

0.98 

2.01 

2.19 

2.13 

2.45 

T.S.  (4.2) 

1.25 

0.99 

1.88 

2.14 

2.30 

2.18 

V 

T.S.  (5.1) 

1.00 

0.99 

2.24 

2.06 

2.44 

2.38 

T.S.  (5.2) 

0.93 

1.05 

1.96 

2.15 

2.20 

2.11 

T.S.  (6.1) 

0.98 

0.84 

1.93 

2.20 

2.34 

2.57 

T.S.  (6.2) 

1.03 

0.90 

2.09 

2.17 

2.16 

2.42 
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ther  investigate  the  relationship  between  ram  and  the  reference 
command  and  OTPV,  the  methodology  of  variance  analysis  is 
employed. 

As  the  test  conditions  are  almost  same,  the  two  results  in  each 
test  state  are  proximately  belong  to  a  Gaussian  distribution.  For 
convenience  of  analysis,  OTPV  is  defined  as  the  factor  A,  and 
the  six  OTPVs  are  defined  as  six  levels,  A\,  Ai,  A3,  A4,  -4.5 
and  A6.  Similarly,  the  mean  value  of  the  reference  command 
voltage  is  defined  as  factor  B.  Because  the  open  voltage  of  the 
gas  flow  regulator  is  about  1.9  V,  the  for  the  voltages  of  1.0  and 
1.5  V  are  meaningless  in  analysis.  Thus,  the  four  valid  results 
for  the  mean  voltages  of  2.0,  2.5,  3.0  and  3.5  V  are  divided  into 
four  levels:  £>1 ,  £>2,  S3  and  S4. 

When  factors  A  and  S  are  at  levels  of  Ai  and  Bj,  Tgfr  can  be 
regarded  as  tgfr,  and  Tgfr  ~  Afijifi  ,  rx2gfr).  In  the  following,  Tgfr 
is  written  as  Xij  for  convenience. 

Let 


SS.4  =  bmJ2(XL.  -  X)2 


SSB  =  gfr YliXj.  -  X)2 

j=l 

a  b 

SSAxB  =  ~  *  j  +  X? 


(17) 


sse=EED%-%)2 


Xij.  =  -J2xijk  i  =  1,  2,  •  ■  ■ ,  6;  j  =  1,  2,  3, 4 

m  k=  1 

;=1,2.3.4 

SJ  1=1  k= 1  1=1 

=ixx> 


In  Eq.  (17),  SS.4,  SSg  and  SS^xfi  are  defined  as  the  sum  of 
squares  for  A,  B  and  Ax  B,  respectively.  In  addition,  SSe  is 
defined  as  the  sum  of  squares  of  the  errors;  a  and  b  are  levels  of 
A  and  S,  respectively;  m,  and  Ax  B  are  sum  of  every  test 
state,  summation  operator  and  interaction  of  A  and  S  on  Xjj, 
respectively. 

Let  n  =  abm. 

According  to  the  theory  of  mathematical  statistics  [27], 
because  X^  ~  J\f(jix,  n2x),  the  statistical  variables,  SS.4,  SSg, 
SS_4xB  and  SSe,  are  independent  of  one  another,  and  the 


statistical  variable  (SSe/tr^)  has  a  y2(n  —  ab)  distribution.  Fur¬ 
thermore,  if  factors  A,  B  and  Ax  B  have  no  effect  on  the 
test  data  Xy*,  the  statistical  variables,  (SSA/crx),  (SSg/cr^) 
and  (SSAxB/:j2x)  have  a  y2{a  —  1)  distribution,  x2(b  —  1) 
distribution  and  y2((a  —  1  )(b  —  1))  distribution,  respectively. 
Furthermore,  (SS^/SSe),  (SSg/SSe)  and  (SS^xg/SSe)  have 
a  T(a  —  1,  n  —  ab)  distribution,  Tib  —  1,  n  —  ab)  distribution 
and  Tiia  —  l)(b  —  1),  n  —  ab)  distribution,  respectively. 

By  analyzing  the  data  in  Table  2  with  the  methodology  of 
variance  analysis,  we  obtain  the  result  listed  in  Table  3,  and 
conclusions  as  follows: 

According  to  the  definition  of  the  ^distribution,  the  critical 
theory  values  of  ^o.95(5,  24),  ^0.95 (3,  24)  and  ^o.95(15,  25)  are 
2.62,  3.01  and  2.11,  respectively. 

The  values  of  T  in  Table  3  are  greater  than  the  critical  theory 
values  of  2^0.95 (5 ,  24),  ^b.95(3,  24)  and  ^0.95(15,  25).  There¬ 
fore,  the  assumption  that  factors  A,  A  and  Ax  B  have  no  effect 
on  Xij  turns  out  to  be  false.  This  result  explains  the  time- varying 
and  nonlinear  property  of  rgfr.  Furthermore,  by  comparing  dif¬ 
ferences  between  the  three  values  of  T  in  Table  3  and  their 
corresponding  critical  theory  values,  ^o.95(5,  24),  ^o.95(3,  24) 
and  7-0.95(15,  25),  it  is  obvious  the  difference  of  factor  A  with 
its  critical  theory  value  stands  out.  This  is  because  factor  A  has 
a  greater  effect  on  rgfr. 


3.  Adaptive  control  strategy  of  airflow 

To  optimize  the  operating  performance  of  the  air-supply 
subsystem,  adaptive  control  is  employed  to  alleviate  the 
dynamic  variation  arising  from  the  time-varying  property  at 
different  operating  points.  This  strategy,  illustrated  in  Fig.  7, 
includes  on-line  parameter  identification  and  real-time  assign¬ 
ment  of  the  closed-loop  poles  according  to  the  identified 
results. 

In  the  process  of  operation,  four  goals  are  to  be  achieved:  (1) 
to  quantify  the  desired  airflow;  (2)  to  define  the  ideal  kinetics 
of  air  replenishment  into  the  cathode;  (3)  to  design  a  con¬ 
troller,  under  the  regulation  of  which,  airflow  responds  to  the 
reference  command  dynamically  and  accurately;  (4)  to  esti¬ 
mate  time- varying  parameters  in  real  time,  modify  dynamics  of 
the  air-supply  closed-loop  subsystem  and  consequently  ensure 
the  airflow  approaches  the  desired  airflow  with  a  certain  preci¬ 
sion. 


71 


Gas  flow 
regulator 


Parameter 

identification  and  pole 
assignment  algorithm 


Fig.  7.  Adaptive  controller. 


Table  3 

Variance  analysis  result  of  rgfr 
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Source  of  variance  Sum  of  squares  Degree  of  freedom  Sum  of  mean  squares  Value  of  F 

A  7.4914  5  1.4983  123.5155 

B  2.8194  3  0.1226  10.1056 

A  x  B  0.5957  15  0.0397  3.2738 

Error  0.2911  24  0 


3.1.  Airflow  supply 


3.3.  Controller  design 


As  airflow  is  directly  related  to  the  output  power  of  a  fuel 
cell  system,  much  attention  has  been  devoted  to  airflow  replen¬ 
ishment.  Gelfi  et  al.  [28]  compared  dynamic  flow  capabilities 
of  a  compressor  of  a  fuel  cell  system  with  a  blower.  The  paper 
showed  the  two  systems  were  dynamically  similar  in  providing 
airflow  in  the  cathode  channel.  The  blower  spent  time  spinning 
its  rotor  inertia,  whereas  the  compressor  spent  time  pushing  air 
and  elevating  supply  manifold  pressure.  Suh  and  Stefanopoulou 
[29]  surveyed  the  dynamic  coupling  between  a  compressor  and 
a  fuel  cell  system  when  the  compressor  motor  was  driven  by  the 
stack  power. 

For  several  reasons  [9,30],  the  airflow  supplied  to  the  cath¬ 
ode  should  exceed  the  airflow  necessary  for  reaction.  This  leads 
to  an  oxygen  excess  ratio,  Aq2,  which  is  defined  as  the  ratio  of 
oxygen  supplied  to  oxygen  used  in  the  cathode.  A  high  oxygen 
excess  ratio,  and  thus  a  high  oxygen  partial  pressure,  improves 
the  output  power  of  a  fuel  cell  stack.  However,  after  an  optimum 
value  of  the  oxygen  excess  ratio  is  reached,  any  further  increase 
will  cause  an  excessive  increase  in  compressor  power  and  thus 
deteriorate  the  system’ s  net  power.  Pukruspan  [  1 2] ,  in  his  doctor¬ 
ate  thesis,  showed  by  experiment  that  the  highest  net  power  was 
achieved  at  an  oxygen  excess  ratio  between  2  and  2.4  depending 
on  the  stack  current  in  the  fuel  cell  system,  whereas  experimen¬ 
tal  data  [31]  showed  the  optimum  oxygen  excess  ratio  of  our 
fuel  cell  system  is  =  1.5.  Therefore  the  expected  airflow  in 
this  paper  is 


/air,ca,ex  =  *02 


Nist 

0.21  x  4 F 


(18) 


3.2.  Desired  kinetics  of  the  air-supply  subsystem 


After  considering  a  variety  of  factors,  the  optimum  dynamics 
of  our  air-supply  subsystem  is 

Tgfr,m  =  —  /air, gfr  +  kwfKgftUgit 


=  /air,ca,ex  ^ 

kwfK gfr 

where  rgfr,m  is  the  optimum  time  constant  of  the  airflow  meter. 
The  corresponding  discrete-time  model  of  Eq.  (19)  is 


/air,gfr(k)  —  ^gfr.n 


i  +  q  1 


1  +  agfrmq~ 


TTjr/air,gff,ex(^) 


where  bgfrm  —  xsj (zs  +  2zgfrm ),  flgfr  —  (zs  2zgfrm  )/(fs  3~ 

2Tgfr ,m) 


To  facilitate  the  design  of  digital  controller,  the  model  of 
the  gas  flow  regulator  Eq.  (9)  is  represented  as  a  discrete-time 
equation  of  the  difference  operator  as 

(1  +  Oair*?- 1 )  /air, gfr  =  &airO  +  q~l)kwfKgtUgflk) 

+fcair(l  +  q~l)e(k )  (21) 


where  bgtl  =  zs/(zs  +  2 rgfr),  flgfr  =  (xs  -  2zgb)/(rs  +  2rgfr); 
and  zs  is  the  sample  time. 

To  ensure  the  airflow  dynamically  tracks  the  reference  com¬ 
mand  with  a  certain  precision  at  a  particular  operating  point,  a 
control  strategy  is  employed  to  assign  poles  of  the  air-supply 
subsystem  for  which  the  time-varying  parameter  is  identified 
on-line,3  and  eventually,  the  dynamics  of  the  air-supply  subsys¬ 
tem  approach  the  described  by  Eq.  (19).  To  achieve  this  goal, 
the  following  control  method  is  adopted: 

W%h.(q-])kwfKg{rUg{r(k)  =  //gfr(*r')/air,gfr;exW 

—  Ggfriq  1 )  /air.gfr  (k)  (22) 

where,  %r(*r')  =  (bgh,m/bam).  Wam{q~l)  and  Gam(q-')  can 
be  solved  according  to  the  following  Diophantine  equation: 

AgfM~l)Wg{M~l)  +  %r(*?-1)Ggfr(<?-1)  =  AgfTim(q-])  (23) 


where 


Agflq  ')  =  1  +  flgfr*?  1 

Bgfriq  1 )  =  bgfr  +  bgftq  1 
Ggfriq  1 )  =  cgff  +  Cgfrq  1 
Agft.miq-')  =  1  +  flgfr.m*?-' 

13gfr,m(q  1 )  =  Hgfriq  1  )Bgfr(q  1 ) 

Wg fr(<?  1 )  =  u>  gfr 
Ggfr(<?  J)  =  ggfr 

Solving  the  Diophantine  equation,  we  get 


WgfrOT1 


tj  _  1  ~  ag 

1  flgfr 


1  Ggff(^  3)  = 


flgfr, m  ~  flgfr 
bgfr(  1  —  flgfr) 


(24) 


3  The  method  for  closed-loop  parameter  identification  will  be  presented  in  the 
following  subsection. 
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Substituting  Wgh(q  '),  Ggb(q  ')  and  Hgu(q  ')  into  Eq. 
(22),  the  control  command  can  be  obtained  as 

UerAk)  =  kwfKgftbait(l  -  a^m)[b™'m{l  ~  flair)/air’gfr’ex(fc) 

—  (flair.m  —  flair)/air,gfr(&)]  (25) 


Furthermore,  the  closed-loop  dynamics  of  the  air-supply  sub¬ 
system  could  be  deduced  in  joint  consideration  of  Eq.  (21)  and 
the  control  command  as 


/air.caOO  = 


W(  1  +  g"1) 
1  + 


/air.gfr.exW 


+ba 


1  ^air.m  1  "h  flair q  '  ,ru\ 
~ j  1  j  _  1 

1  —  flair  1  +  «air,m<7 


(26) 


then  Eq.  (27)  can  be  reformed  as 

Pgfr  =  WgfrCg(Mk)  (30) 

where 

{■Pgfr  =  Pgfr.O  +  Pgf r,l?_1 

Pgfr.O  =  /air,gfr(&)  —  ^gfr/air.gfr.ex^)  (31) 

Pgfr.l  =  agfr/air,gfr(&)  -  %r/air,gfr,ex(fc) 

According  to  the  principle  of  parameter  identification,  if  Pgu 
and  Cgfr  are  prime  to  each  other,  elements  of  Pgfr,  pgh,o  and 
pgfr,i  are  identifiable  as  e(k')  belongs  to  Gaussian  distribution. 
As  the  only  factor  of  Cgfr  is  (1  +  <y~'  ),  to  prove  Pam  and  Cam  are 
prime  to  each  other  is  to  prove  there  is  no  (1  +  q~l)  in  Pgfr.  As 
deg  Agb{q~l)  =  1,  if  pghfi  /  Pam,  l,  it  can  be  concluded  there 
is  no  (1  +  q~l )  in  Pam. 

Because 


3.4.  Closed-loop  parameter  identification 


Pgfr.l  —  Pgfr.O  —  (flgfr  —  l)/air,gfr(^) 


(32) 


Variation  of  the  main  time  constant  of  the  gas  flow  regula¬ 
tor  mainly  arising  from  its  time- varying  property,  often  occurs 
with  the  shift  of  the  operating  point  in  the  air-supply  subsystem. 
This  directly  leads  to  a  discrepancy  between  the  real  airflow 
dynamics  and  its  desired  counterpart.  To  reduce  the  influence 
of  the  variation,  on-line  modification  of  dynamics  of  the  air- 
supply  subsystem  is  necessary.  The  premise  this  application  is 
the  identification  of  the  main  time  constant  in  the  closed-loop 
air-supply  subsystem.  As  is  well  known,  parameter  identifica¬ 
tion  in  a  closed-loop  system  is  more  difficult  than  that  in  an 
open-loop  system,  and  in  most  cases,  this  task  is  unrealizable. 
Next,  we  will  prove,  in  theory,  that  time- varying  parameters  in  an 
gas  flow  regulator  are  identifiable  in  the  closed-loop  air-supply 
subsystem. 

Proof.  Considering  the  model  of  the  air-supply,  i.e.  Eq.  (21) 
and  its  control  equation  Eq.  (22),  the  closed-loop  subsystem  can 
be  expressed  as 

(Agfr{q-])Wgh(q-')  +  %r(4“1)Ggfr(tf“1))/air,gfr(A0 
=  Hgfr(q- 1  )Bg\\(q- 1  )/air,gfr,ex(k)  +  CgfT(q-])WgfT(q-]Mk) 

(27) 


it  is  easy  to  conclude  that,  if  agfr  A  1 ,  then  pgu/)  ^  Pgfr.i  •  Due  to 
agfr  =  (ts  —  2rgfr)/(rs.  +  2rgfr),  it  is  obvious  that  the  inequality 
of  agfr  /  1  comes  into  existence,  namely,  pam, o  /  Pam,  l  ■  It  can 
be  deduced  the  elements  Pam  and  Cam  are  prime  to  each  other. 
So  far,  conclusion  can  be  deduced  that  the  elements  of  Pam, 
namely  pam, o  and  pam, l  are  identifiable. 

Assuming  that  a  certain  method,  such  as  least  squares  identi¬ 
fication  or  maximum  likelihood  estimation,  is  applied  to  identify 
pgfro  and  Pgfr, i .  the  results  are  pg fr>o  and  pgh,i,  respectively. 

In  addition,  because  bg fr  =  rs/(2ram  +  zs)  and  aam  =  {ts  — 
2 tam)/{ts  +  2ram),  it  is  clear  that  bam  and  aam  are  dependent  on 
each  other,  and  their  relationship  can  be  expressed  as 

2Agfr  =  agfr  +  1  (33) 

Applying  pgbfi,  pg fti,  Hgb(q~l)  =  (bgfgm/bgfr),  2 bam  = 
aam  +  1  and  substituting  Eq.  (24)  into  Eq.  (29),  we  can  get 

(  Pgfr.O  —  /air,gfr(^)(l  +  agfr,m )  N 
+^gfr,m/air,gfr,ex(^)  ^4) 

Pgfr.l  —  2agfrm /air  gfr(k) 

\+Vr,m/air,ca,exW  / 

where 


( ^gfr,//air,gfr,ex (k)  2fa\r  g(r(k)  +  Pgfrfl  (flgfr.m  l)/air,gfr(£)  \ 

y  ^gfr.m /air.ca.ex  (k)  —  /air.gfr(^)  —  Ugfr.m /air.gfr(^)  +  Pgfr.l  (agfr,m  —  l)/air,gfr(^)  J 


where 

deg  Wgft(q  ])  =  0 

deg  Ggfr((7-1)  =  0  (28) 

dLe|Wgfr  (</-')  =  0 

Pgfr  =  —  Hgfriq  )Pg(Aq  )/air,gfr,ex(£)  +  ( PgiAq  )Ggfr(q  ) 

+  Agb{q~l)WgfM~X))f™,&{k)  (29) 


and 

det(M)  =  Pgfr.l  -  Pgfr.O  -  (agfr,m  -  l)/air,gfr(&)  (35) 

where  det(M)  is  the  determinant  of  the  square  matrix  M. 

According  to  Eq.  (32),  if  the  dynamics  of  the  air-supply  sub¬ 
system  converge  to  Eq.  (20),  the  determinant  of  M  converges  to 
zero.  This  means  the  two  equations  in  Eq.  (34)  are  linearly  depen¬ 
dent  on  one  another,  and  the  estimation  of  agfr  can  be  obtained 
by  solving  either  of  them.  At  the  same  time,  an  estimation  of 
bgf{  can  be  inferred  from  Eq.  (33). 
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Next,  least  squares  identification  is  discussed  to  estimate  agfr 
and  figfr  in  the  air-supply  closed-loop  system. 

It  is  easy  to  rearrange  Eq.  (21)  into  an  ideal  form  for  least 
squares  identification  as 

/air,gfr(fc)  =  hT6  +  bgfrwgfr(e(k)  +  e(k  -  1))  (36) 


H DC/DC  \~ 


CuITCnt  p|  Load  I 

regulator  «— i - 1 - 1 


where 

f  h  =  [0.5kwfKgft(u^k)  +  ugfT(k  -  1)) 

\  — /air,gfr(^  —  1),  0.5kwf  KgfAUghik)  +  Ugfr(k  —  1))]^ 

(  e  =  [agfr.  If 

In  addition,  in  the  air-supply  subsystem,  /air.gfr  can  be  mea¬ 
sured  directly  by  the  flow  meter  built  in  gas  flow  regulator. 

A  recursion  algorithm  can  be  expressed  as 

r  Egfr(O)  =  0!gfr/2x2 

9(0)  =  pg  fAk  +  1)  =  Kgfr(k  +  i)\uo(k  +  1) 

+m)-hTgfI(k+mk)] 

<  Kg(r(k  +  1)  =  [hTg(r(k  +  \)Pgfrhgfr(k  +  1)  (37) 

+lLgfrr'  Pgfr(k)hgfr(k  +  1) 

PgfAk  +  1  )=-[/-  Kgir(k  +  \)hTg{r(k  +  V)]Pg(Ak) 

P-gfr 

0  <  /Tgfr  <  1 

where  agfr  is  a  sufficiently  large  real  number,  and  Pg(T  is  a  suffi¬ 
ciently  small  2x1  dimensional  vector. 

According  to  Eq.  (33),  the  estimation  of  bg fr  could  be  deduced 
as 

,38, 

In  the  air-supply  subsystem,  it  is  practical  to  estimate  the 
parameters  of  agfr  and  bgfr  by  least  squares  identification  and 
subsequently  to  modify  the  kinetics  of  the  closed-loop  air-supply 
subsystem.  This  method  is  the  so-called  adaptive  strategy  for  air 
replenishment  control  of  the  air-supply  subsystem.  The  step-by- 
step  procedure  for  this  strategy  is  expressed: 

(1)  Initialize  figfr(0). 

(2)  Get  uigfr,  ggfr  and  hgfr  by  solving  Eq.  (23). 

(3)  Determine  /air.gfr  by  solving  Eq.  (22)  and  subsequently 
apply  it  to  the  air-supply  system. 

(4)  Save  /air,gfr(k  —  1)  and  sample  /air.gfr(^),  and  at  the  same 
time  determine  the  reference  command  ugfr(k). 

(5)  With  /air,gfr(^  1),  ,/air, gfr(^)  and  Ugf, ,  estimate  cigf^(k)  and 

bgfr(k)  by  least  squares  identification. 

(6)  Return  to  step  (2). 

4.  Test  and  experiment 

The  fuel  cell  test  bench  in  this  paper  was  designed  and  assem¬ 
bled  by  CASIA  (Institute  of  Automation,  Chinese  Academy  of 
Sciences).  This  test  bench,  as  shown  in  Fig.  8,  mainly  includes 
the  fuel  cell  system  described  above,  a  dc/dc  converter  for  stack 
voltage  stabilization  (dc/dc  output  voltage  of  12  V),  a  current 


Fig.  8.  Scheme  of  test  bench. 

regulator  for  stack  current  control,  a  battery  as  auxiliary  power, 
a  regulatable  resistor  as  a  load  with  maximum  power  of  800  W 
and  a  electric  motorcycle  with  maximum  power  of  500  W.  To 
verify  the  control  strategy  described  above,  a  test  comparison  is 
made  under  the  control  strategy  of  the  pole  assignment  described 
in  Section  3.3  with  and  without  rgfr  identification.  In  the  former 
case,  the  ideal  value  of  rgfr  is  set  at  2.3  s,  and  in  the  latter  case, 
the  value  of  rgrr  varies  according  to  the  shift  of  the  operating 
point. 

In  tests,  operating  points  of  the  fuel  cell  system  are  regulated 
through  changes  of  airflow  in  the  gas  flow  regulator  and  of  stack 
current. 

(1)  Test  1 

The  stack  current  is  regulated  to  reach  40  A  by  coor¬ 
dination  of  the  current  regulator  and  the  resistor,  and  the 
cathode  pressure  is  controlled  at  30  kPa  (relative  to  atmo¬ 
spheric  pressure).  There  are  two  airflow  operating  points  at 
0.00176  and  0.00238  mol  s-1.  Step  disturbances  of  airflow 
at  0.0002  mol  s_1  are  exerted  on  the  gas  flow  regulator  as 
reference  command  changes  at  these  two  operating  points. 

Test  data  are  shown  in  Figs.  9  and  10.  Theoretically,  the 
time  span  for  the  step  response  of  a  first  order  linear  system 
to  reach  61.9%  of  the  stable  output  altitude  is  the  system’s 
time  constant.  The  dotted  lines  in  Figs.  9  and  10  represent 
61.9%  of  the  stable  output  airflow,  and  thus,  the  abscissa 
of  the  intersection  point  between  the  dotted  line  and  test 
data  (solid  line)  is  equal  to  the  time  constant  of  the  gas  flow 
regulator. 


Fig.  9.  Test  data  at  the  operating  point  of  /an.gfr  =  0.0176  mol  s  '  . 
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Step  response  of  airflow  without  parameter  identification 


Step  response  of  airflow  without  parameter  identification 


8  9 

Time(s) 

Step  response  of  airflow  with  parameter  identification 


Time(s) 

Fig.  10.  Test  data  at  the  operating  point  of  /ari,gfr  =  0.0238  mol  s“ 


2  3  4  5  6 

Time(s) 

Step  response  of  airflow  with  parameter  identification 


Fig.  12.  Test  data  at  the  operating  point  of  ist  =  40  A. 


It  is  easy  to  see  the  time  constants  in  the  system  with 
parameter  identification  at  the  two  operating  points  are 
almost  equal  at  2.3  s.  However,  the  time  constants  in  the  sys¬ 
tem  without  parameter  identification  vary  obviously  from 
1.5  to  2.4  s. 

(2)  Test  2 

The  OTPV  is  fixed  at  a  suitable  position  and  the  stack 
current  operating  points  are  fixed  at  30  and  40  A  with  cor¬ 
responding  airflows  of  0.00133  and  0.00237  mol  s-1.  Step 
disturbance  of  the  stack  current  is  5  A  at  these  operating 
points.  The  controller  can  adjust  airflow  to  track  the  step 
disturbance  of  the  stack  current  with  an  oxygen  excess  ratio 
of  1.5.  The  dynamics  of  the  gas  flow  regulator  are  compared 
for  these  closed-loop  air-supply  subsystems. 

According  to  the  test  data  in  Figs.  1 1  and  12,  it  is  obvious 
the  time  constants  in  the  system  with  parameter  identifica¬ 
tion  at  the  two  operating  points  are  almost  equal  at  2.3  s. 
However,  time  constants  in  the  system  without  parameter 
identification  vary  obviously  from  1.77  to  2.1  s. 


Step  response  of  airflow  without  parameter  identification 


flow  with  parameter  identification 


— 

Fig.  11.  Test  data  at  the  operating  point  of  ist  =  30  A. 


(3)  Test  3 

The  motorcycle  is  connected  to  the  current  regulator 
as  load  and  the  initial  power  of  the  cycle  is  controlled  at 
about  120  W  which  responds  to  stack  current  at  about  9  A. 
Abrupt  acceleration  is  performed  through  regulation  of  the 
accelerograph.  Stack  current  change  to  15.5  A  in  the  adap¬ 
tive  controller  case  and  15  A  in  the  non-adaptive  controller 
case,  respectively.  Test  data  in  Fig.  13  shows  the  time  con¬ 
stant  is  about  2.3  s,  i.e.  the  ideal  value,  in  the  adaptive 
controller  case,  while  the  time  constant  is  about  1.8  s  in 
non-adaptive  controller  case  which  obviously  departs  from 
the  ideal  one.  From  the  test  results  in  Figs.  9-13,  it  can 
be  concluded  adoption  of  the  adaptive  control  strategy  can 
be  used  to  adjust  the  gas  flow  regulator  time  constant  to  a 
fixed  point,  and  consequently  make  the  gas  flow  regulator 
perform  with  an  desired  dynamic  action. 


Time(s) 


Fig.  13.  Test  data  in  motorcycle  i 
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5.  Conclusion 

Technically,  PEM  fuel  cell  systems  could  become  viable 
propulsion  system  alternatives  to  conventional  combustion 
engines,  as  long  as  a  sufficient  amount  of  power  output  is  avail¬ 
able  to  meet  the  need  of  the  load.  To  achieve  this  goal,  every 
subsystem  of  a  fuel  cell  system  has  to  act  smoothly  and  be 
dynamically  co-operative  with  each  other.  However,  the  time- 
varying  and  nonlinear  characteristics  of  these  subsystems  are 
obstacles  during  operation.  To  run  the  air-supply  subsystem  to 
meet  the  need  of  air  replenishment  into  the  cathode,  this  paper 
proposes  an  adaptive  control  strategy  to  regulate  the  dynamics 
of  the  subsystem  at  different  operating  points.  A  summary  of 
new  findings  obtained  is  as 

(1)  Experimentally  find  the  time- varying  and  nonlinear  proper¬ 
ties  of  the  gas  flow  regulator. 

(2)  Optimize  the  kinetics  of  the  air-supply  subsystem  by  design¬ 
ing  an  algorithm  for  pole  reassignment  via  the  solution  of  a 
Diophantine  equation. 

(3)  Theoretically  find  the  the  time-varying  parameter  of  the  gas 
flow  regulator  identifiable  in  the  closed-loop  system,  and 
subsequently  apply  the  least  squares  identification  algorithm 
to  estimate  the  time- varying  parameter  and  eventually  pro¬ 
pose  an  adaptive  strategy  for  the  control  of  the  air-supply 
subsystem. 

(4)  Compare  the  control  strategy  with  and  without  parameter 
identification,  and  verify  the  validity  of  the  adaptive  strategy. 

The  tests  have  shown  the  designed  adaptive  control  strategy 
can  lead  to  good  system  performance  and  stability.  Although 
this  paper  is  limited  to  the  control  of  the  air-supply  subsystem 
only,  it  is  important  to  note  the  adaptive  principles  can  be  easily 
applied  to  other  subsystems  of  the  fuel  cell  system. 

Further  work  will  include:  (1)  examination  of  the  influence  of 
cathode  pressure  on  airflow,  (2)  investigation  of  the  stack  current 
impact  on  the  kinetics  of  air  replenishment  into  the  cathode  and 
(3)  integration  of  an  air-supply  subsystem  with  other  subsystems 
into  an  organic  unity. 
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